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^ ■ Abstract 

QO ■ A c-number path integral representation is constructed for the solution of the 

Dirac equation. The integration is over the real trajectories in the continuous 
\ three-space and other two canonical pairs of compact variables controlling the 

i> ■ spin and the chirality flips. 

The path integral representation of the quantum amplitudes has become an essential 
device in the development of an intuitive picture about quantum physics, the organisation 
of the perturbation expansion, and in establishing non-perturbative approximations. A 
serious limitation of the functional methods arises from the difficulties in extending the 
path integral formalism for relativistic fermions. There have been numerous proposals to 
solve this problem. The original attempt at a path integral representation for the solution 
of the 1+1 dimensional massless Dirac equation was a random walk along the line cone of 
a space-time lattice with a Poisson distributed helicity flips The method was soon 

generalised for 3+1 dimensions M and the highly complicated time ordered products of the 
Dirac matrices was dealt with by the integration over Grassman- valued trajectories J3]. The 
^ propagation along the space-time trajectories and in the spinor space was separated and the 

latter was treated by brute force in refs. 0. The propagation can be described in a five 
dimensional manifold and the resulting path integrals resemble the proper time formalism 
||. Guided by supersymmetry J7[ the Grassman variables were introduced combined with 
the proper time ||. We aim at a path integral formalism which stays as close as possible 
to the physics of the problem (does not involve proper time, the continuous trajectories are 
parametrized by the real time), and is based on c-number integration to allow the use of non- 
perturbative methods, such as the semiclassical approximation for the fermion propagator. 
Such a formalism, employed with the second quantization might bring bosons and fermions 
on equal footing. 

The goal of the present paper is the path integral solution of the first quantized Dirac 
equation 

ihd t ^ = H D m = [(-ichd + eA) ■ a + mc 2 (3 - eV] (1) 

where *&(x, t) denotes the Dirac bi-spinor, (A(x,t),V(x,t)) is an external gauge field and 
a i — 7o7j; P — 7o are Dirac matrices. Before embarking the technical problems to resolve 
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in finding the path integral one is faced with the usual inconsistency of the first quantized 
quantum mechanics, namely the existence of states with negative energy. This appears in our 
context as the ill defined nature of the time evolution operator g-^^W^ even for imaginary 
time. We avoid this problem by the replacement H D — > H cuto ff = H D + p 2 /2M in the 
exponential. The mass parameter M introduced in this manner plays the role of an UV 
cutoff because this modification is acceptable only for p < Mc. The cutoff is implemented 
in the path integral formalism by the small but finite time step At. We shall find the 
relation Mc 2 = nh/At where k is a real positive constant, indicating that the time evolution 
operator is well defined when Im(tAt) < 0, for the time axis which is tilted as t — > te ta , 
— 7i / 2 < a < 0. The choice a = —it /A will serve as the safe starting point to establish 
the convergence of the path integral expressions introduced below by means of the analytic 
continuation either to real or imaginary time. 

The problems to solve in the path integral formalism is (i) to generate the first order 
derivatives in the space coordinates and (ii) to treat the spinor indices. The problem (i) is 
circumvented within the framework of the Wiener-type path integration by using the bare 
cutoff hamiltonian, H cuto ff- The main difference with the non-relativistic case is the unusual 
cutoff dependence || . (ii) is solved by a well known method, the introduction of additional 
continuous degrees of freedom for the representation of the dynamics in the Dirac space [1TDI . 



In the first step we construct the integration over the three-space trajectories by keeping 
the spinor dynamics in its algebraic form. For this end we introduce the infinitesimal time 
propagator G(x, y; At) = e lA ^ x ^ T ) defined as 

V(x, r + At) = J dyG(x, y; Ar)*(y, r), (2) 

where r = tc. One may motivate the choice of A(x, y; At) for free massless fermions in the 
following manner. The self reproducibility, 



e iA(x,y,T 1 +T 2 ) 



J dze iA{x ' z ' n) e iA ^' mT2 \ (3) 



can easiest be satisfied by using quadratic expression in the coordinates for A(x, y; At). The 
translation and rotation symmetries require the functional form 



A(A 



T 



A(x, y, At) = ^p(* - yf + B(At)(x - y ).S + C, (4) 

Sj standing for a vector which transforms as the spin under rotations. Since the equation of 
motion contains no U we keep A fr- independent. Then the choice A = a/ At 2 and B = h/ At 
follows by dimensional argument. We can carry out the integral in (0) after expanding the 
amplitude in the right hand side at y = x one arrives at 



*(i,r + Ar) = J dye l ^ {x ~ y)2 -^ {x - y) ^ +c h{y,t) 



a 



tt(x,t) (5) 



for the appropriate choice of the constant C and having dropped the contributions 0(At 2 ). 
This agrees with the massless Dirac equation for Sj = ctj, a = b = k if the characteristic 
length of the amplitude is long enough, At << fr/\p\, to ignore the terms 0(At 2 ). The 
mass represents a mixing between the chiral spinors and can be taken into account by the 
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modification G — ► Ge 1 / 3At / a ; A = Ti/mc and Ar << A. The resulting path integral in the 
absence of the external field is 



2 l At J At 



+iC*(Ar) 



-iArS£/3 



(6) 



for a suitable chosen C(At). One introduces the formal notation lim TL f dx«e lC ^ r ^ 



J and can writes 

G?(x,y;r) = |l?[x]r 



1_ /" dx \ dx 

2 \dr ) dr' 



(7) 



The convergence of the path integral can be established by starting with the time axis tilted 
by — 7r/4 as explained in the introduction and performing the analytical continuation for 
real time. The path integral can easily be computed for the massless case, 



G(x,y;r) 



d 3 p 



At j 

2k 1 



-I- 



+c 



(8) 



in agreement with the form of H cuto ff. 

The typical trajectory in (0) is more regular than in the non-relativistic path integral 
because the linearly diverging coupling constant At -1 of the action suppresses the fluc- 
tuations. The divergence in the exponent of the integrand is not a real surprise, one is 
accustomed to find diverging parameters in the bare, cutoff quantum field theories and even 
in non-relativistic quantum mechanics ITT]] . In fact, the cutoff of a renormalizable model is 



removed in the following manner: We are given first a regulated path integral involving the 
lagrangian L(4>, <9 M 0; g n ) and a cutoff A. The cutoff is then sent to the infinity. The parame- 
ters of the lagrangian become functions of the cutoff, g n —>■ g n (A) in such a manner that the 
regulated path integral with the lagrangian L(<f>, d^; g n (A)) yields convergent observables 
as A — >• oo. The fine tuning of the bare parameters is introduced in our case for the same 
goal, to arrive at a cutoff independent, convergent result, (0). One can verify the scaling 
relations 



-V) 

At) ' 



{a/At 



Ar << A, 
Ar >> A 



(9) 



confirming the relativistic crossover at Ar m A . These scaling laws indicate that the non- 
differentiable fractal trajectories of the non-relativistic quantum mechanics are "regulated" 
by embedding them into a smoother relativistic propagation. 

In the second step we eliminate the time ordered product of matrices in the spinor space 
UJ . For this end we write the matrices of the Dirac hamiltonian in the chiral representation 



as the direct products ctj = uj ® r 3 and (3 — Ti, where a and r are the Pauli matrices acting 
in the spin and the chirality spaces, respectively. In the matrix formalism (^]) one uses the 
basis vectors \m— |) m = 0, 1 in the spin and the chirality space. Instead of these vectors the 
dynamics in the four dimensional Dirac space will be followed by means of the overcomplete 
and continuous basises |0) and \x), introduced for the spin and the chirality, respectively, 
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«> = 4= E e— |m-i), (10) 



m=0,l 

with — 7r < a < 7r. The bi-spinor ty(x,t) is converted into a single component amplitude by 
projecting it onto the basis vector \<f>,x) — \4>) ® \x)-> ^{4>i X; x , t) = (0, x\^( x , t)). The time 
evolution is given by 

X, a;, t + At) = f dyd^'dx'G(x, 0, x, y, <fr', X\ Ar)tf x', 1/, r), (11) 

where 

G(x,0, X ,l/,0',x'; Ar) = (<f>, X \e-i ATH ^\<l>',x')- (12) 
By the help of the relation 

7T 

f 1 1 

y d0(m - -|0}(0|m - -) = 5 m /, m , (13) 
m, m' = 0, 1 we find the resolution of the identity 

da | a) (a | = 1. (14) 



The repeated insertion of this form of the identity into the time evolution operator yields 
the path integral where the trajectories are x(r), <p(r) and x(t). The Poisson resummation 
formula 

E/W = E / dpfipy 2 ™*, (is) 



n=— oo„ 



where n are integers and < e± < 1 can be used to write the overlap between the 

basis elements as 



(a\(3)= f] / g^-^n)^ (lg) 



suggesting the interpretation of the variables a and p as canonical pairs. 
We need in (|5]) the matrix elements of the Pauli matrices, 

gi(a-0) _ i 

(a|(7 3 |/3) = — 

m " J=0,X 



pg*p(a-/3+27rn) 

" ~ ~~ x —I) 



« 7° d P / i 
2 s / yA p ~2 



ip(a— f3+2wn) 



-0 
1-0 



2 V / ^. pe i(p+h)(<*-P+*™)^ (17) 
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where the Poisson resummation formula was used in the third equation to write the sum 
over m as an integral. The matrix elements of a± = o\ ±i&2 are given by (a\a + \f3) = e ia /iT, 
(a\(T-\/3) = e~ ip /n, 



(a\a ± \(3) = -^S nfi e^ +1 ^-^^f(n) 



(18) 



where f(n) is an arbitrary function satisfying /(0) = 1. The Poisson resummation formula 
yields 



1-0 



ft 



n=-oo_ 1+Q 



(19) 



for the choice f(n) = yl — n 2 . The application of ([T7|) and (|1^) gives the matrix elements 
of v ■ a = v 3 a 3 + (u + cr_ + V-<j+)/2, v± = v 1 ± iv 2) 



1-0 



(a \v ■ a\ 



2 E 



-1+0 



^E e »(P+|)(a-/3+27rn) 
27T 



x i v 3 p + \Jl-p 2 



V\ COS 



ft + /3 + 27m 



+ v 2 sin 



'ft + /3 + 27m' 



(20) 



We are now in the position to rederive (|7|) in the continuous basis, 

ds 3 j dt 3jj 



2V-1 , . JV 

G{xf,<f>f,Xf,Xi,<f>i,Xi',r) = II ( / dx k d<f) k dxk) Yl 



fc=l 



2tt 2tt 



E 



(21) 



as iV — > oo, with 



^ = E[( s 3,i + ^) - 0,-1 + 2^) + (t 3>j + ^) ( Xi - Xi-i + Ztt^-) 



+- 



Ar 



-4Ar— t 3J s j • 



I 1 



9 Ar 



(22) 



+ Ar^-(x i )+C 



x = Xi, 0o = (f>i, Xo = Xh %n = x f , (f) N = 0/, Xn = Xf, Aj(x) = A(x,jr/N), Vj(x) = 
V(x, jr/N). We used the mid-point prescription in the vector potential to ensure gauge 
invariance and the vectors corresponding to the spin and the chirality spaces are given by 



s 3 — s 3,j s 



1 _ s 2 §(^+^-1+2^) 
h = tsj U = V /r^|- e 4( % +X,-i + 2 7r iV,)_ 



(23) 



By the shift of the variables <pj — > <pj — 2n J2k<j n k, Xj ~* Xj ~ 2tt J2i<j and the extension 
of the region of integration for cf) and x over the whole real axis the summation over rij and 
Nj decouples and can be ignored. The resulting path integral expression is 
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G ( x f, <Pf, Xf, Xi, <t>u Xi) r ) = II (/ dx k d(j) k dxk^ II ( f e lA 



2tt 2tt 



(24) 



where 



N 
3=1 



S3J + 



+ (*3,j + g J (Xj -Xj-i) 



2 V Ar 



3J Ar 



Ar 



(25) 



-4Ar^-.A,(^±^) + Ar-|^(,,) + ^ 



and 



V 



Sgj COS w - \ 



1 _ s 2 . gin 



, tj 



^/r^|. cos xi±xti 



v 



1 _ £ | . sin ^±1^1 
*3J 



(26) 



/ 



The integration over the variable <fik, Xk can be restricted within the interval [— 27r, 27r] due 
to the periodicity of the integrand. One thus finds in the continuum limit the expression 



D[x)D[s)D[t]e 



(27) 



for the solution of the Dirac equation, where the trajectories s(r) and t(r) are on the sphere 
of unit length. 

Several remarks are in order at this point, (a) As announced above, D[s] = D[cos6]D[<f)] 
and D[t] = D[cosQ]D[x] represent the integrations over the orientation of the spin 
and the chirality which are parametrized as s = (sin 9 cos 4>, sin 9 sin 0, cos 8), and t = 
(sin cos x, sin B sin x, cos 0), respectively, (b) Note that the half of the winding num- 
bers \ J d(f) and | / d\ in the action generate the sign —1 for rotations by 2n. Thus the path 
integration involves double valued amplitudes in the spin and the chirality spaces. This is 
equivalent by saying that the integration D[s] and D[t] is over the covering space which con- 
tains S3 twice, in agreement with the mid-point prescription for s and t, i.e. the apparence 
of the half angles in (p6|) . (c) The variables (cos^, 0) and (cos0,x) can be considered as 
canonical momenta and coordinates. The continuous formal notation actually hides the 
fact that there are one less coordinate than momentum variable in the path integral. The 
terms s^dip/dr and t 3 dx/dr are of the usual formpg, characterizing the path integrals in the 
phase space. The discrete spectrum of the angular momentum and the corresponding finite 
integration range for S3 and t% prevent us from eliminating the momentum variables and 
to reduce the amplitude to a path integral in the coordinate space only. This complication 
leaves behind a rather singular trajectory structure for the angle variables, (d) The term 
involving t\ generates the rotation of the chiral angle 0. In the massless case the \ integra- 
tion gives dt^jdr = 0, the conservation of the chiral angle for massless fermions. The second 
and the fourth terms in the exponent are responsible for the spin precession modulated by 
the chiral angle, (e) Though the parameter k is left free we recall that the saddle point 
trajectory x c i(t) of the path integral (0) is the smoothest for k = tt |[L2|| . (f) It is worthwhile 
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mentioning the double role the piece 0(x 2 ) of the action plays: It generates a second order 
finite difference equation for the saddle point trajectories as expected from the semiclassical 
limit and eliminates the species doubling by appearing as a Wilson term. Note the manifest 
chiral invariance. 

Finally, a remark concerning the second quantization, mentioned at the introduction. 
Due to the spin-statistics theorem the phase corresponding to the rotations by 2tc and the 
exchange of equivalent particles are the same. Since the factor —1 of the rotations is built in 
(p6|) one hopes to arrive at the second quantized path integral without Grassman variables. 
Furthermore, the chiral invariance of the first quantized amplitude (^7|) suggests that the 
no-go theorem for lattice chiral fermions [13] could be avoided. 
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